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Abstract 

Since the discovery of X(3872), its structure has been in ceaseless dispute. The data of X(3872) — )• 
Tr~^7r~Tr'^ J ftp suggest that X(3872) may be a high-spin charmonium-hke of 2 ^. In terms of the 
hght front quark model (LFQM) we calculate the rates of the radiative decays X(3872) — )■ J /ip{ip')^ 



CN I supposing X(3872) to be a 2 ^ charmonium. Within this framework, our theoretical prediction on 

flj \ B1Z{X{?>^72) — 7- il){lS)^) is at order of 10~^ which is slightly lower than the Babar's data but close 

Q . to the Belle's. Our prediction on i37^(X(3872) — ;• -(/''t) is at order of 10"^ if ■0' is a pure 2S state 

CN , or 10~^ if ijj' is a 25 — ID mixture, which does not conflict with the upper bound set by the Belle 

_, collaboration, but is much lower than the Babar's data. Thus if the future measurement decides 

Qh' the branching ratio of BTZ{X{'iS72) — )• ijj'^) to be much larger than 10~^, the 2 '" assignment for 

Ph' X(3872) should be ruled out. 
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I. INTRODUCTION 



Very recently, a series of charmonium-like resonances has successively been observed by 
various experimental collaborations, such as X(3872)[l|, X (3940) [2], F(3940)|3|, and several 
bottomoniun-like states were also discovered at the Belle and Babar energy ranges, such as 
Z(4430)^[4], Zb and ^^ [sl]. It is noted that there is almost no room in the ground state 
representations oiO{'i)®SUf{?>)®SUs{2) to accommodate those newly observed resonances. 
Not only their mass spectra, but also their behaviors of production and decay may hint if 
they are exotic states such as hybrids, molecular states, tetraquarks, or just radial and/or 
orbital excited charmonia. 

It would be an interesting task to determine their inner structures from both experimental 
and theoretical aspects. As a matter of fact, there are different interpretations for these 
resonances which are of different constituent compositions from the regular mesons, so that 
for clarifying its structure, all the interpretations should be studied one by one and then we 
can see if the predictions are in agreement with data. The most reasonable way is to calculate 
its mass and decay widths by assuming its structure in certain theoretical frameworks, then 
a comparison of the prediction with data would confirm or negate the assumption about the 
identity of the resonance. 

Among the newly observed resonances, X(3872) which was found a while ago has caused 
special interests of experimentalists as well as theorists |6l-lll|. Some authors regard it as 
a molecular state 12|-|l4l|. whereas other groups consider it as a tetraquark|15Nl7l| . Instead, 
the recent data on X(3872) -^ tt'^tt^tt^J/iP may hint a possible assignment ^D2 ( i.e. 2~"^) 
for X(3872) which has inspired theoretical interests 2lN24l]. The charmonia mass spectrum 
has been calculated in the pot ential modeled the theoretical prediction for a ^D2 charmo- 
nium is about 3.81GeV 



25l . 1261 ] or 3.84GeV 27| which is 30-60 MeV lower than the measured 
value of X(3872). Because of uncertainty existing in the calculation of the binding energies 
of excited states with the potential model such difference may not be too serious, but the 
difference might imply that our understanding on its assignment might be not completely 
correct. To determine its structure one needs more information from other sources, especially 
its decay modes. The Belle and Babar collaborations also reported their measurements on 



the radiative decays of X(3872)[18|, l28|]. The data given by the Babar collaboration [29| 
are Bn{B^ -^ X (3872) K^)Bn{X (3872) -^ J/ij-f) = (2.8 ± 0.8{stat) ± 0.1{syst)) x 10-^ 
BTZ{B^ -^ X(3872)i^±)S7^(X(3872) -^ tfj'-f) = (9.5 ± 2.7{stat) ± 0.6{syst)) x lO'^ 

3.4 ± 1.4. The Belle collaboration also reported their new re- 
^ X (3872) K^)Bn{X (3872) -^ J/tfj-f) = {1.78l^oil ± 0.12) x 10"^ and 

~^ and 
The 



'^'^^^ Bn(X{3872)^J/ip-y) 

suh as Bn{B^ 



set an upper bound Bn{B^ -> X (3872) K^)Bn{X (3872) -^ ^p'-f) < 3.45 x 10 



B7i(X(3S72)-^M ) ^ ^.1 because they found no evidence for X(3872) — )■ ■;/'S|30[|. 
data from the two collaborations are consistent on X(3872) — t- J/ip'-f but largely apart 
on X(3872) -^ ip'^. Obviously the further measurement and theoretical study are badly 



needed. 

In this work we investigate the radiative decays of X(3872) which is supposed to be a 
^D2 charmonium in the hght-front quark model. The results may help us to determine the 
structure of X (3872). 



The light-front quark model(LFQM) is a relativistic model 3l|, |32| and in this model wave 
functions are manifestly Lorentz invariant and expressed in terms of the fractions of internal 
momenta of the constituents which are independent of the total hadron momentum. This 
approach has been applied to study many processes and thoroughly discussed in literatures 



44j . Generally the results obtained in this framework qualitatively coincide with the 
experimental observation on the concerned processes and while taking the error ranges into 
account (both experimental and theoretical), they can be considered to quantitatively agree 
with the available data. 

To evaluate the transition rate in the LFQM one needs to know the wave functions of 
the parent and daughter hadrons. The wavefunctions for the s-wave and p-wave were given 



in Ref . (37| and we studied the wavefunctions of the d-wave in Ref . J45| with which we are 
able to investigate the transitions involving s-, p- and d-wave mesons. Then for obtaining 
the transition amplitude, one also needs calculate those form factors in terms of the effective 
Lagrangian. In this work we will deduce the form factors for the radiative decay of 2 ^ — )■ 1 
in the covariant light-front quark model. With those formulas we calculate the rate of the 
radiative decays of X(3872) which is supposed to be a pure 2~^ meson. Because of the 
so-called p — tt puzzle, namely the branching ratio of tp' — )■ pir is extremely small while 
J/ip — ^ pvr is one of the main decay modes of J/ip, it is suggested iJj' might not be a pure 2S 



state {i/j{2S)) but a mixture of 25* and ID states {ip{2S — ID) )|46(|. This allegation has not 
been fully confirmed so far, even though it forms a reasonable interpretation for the p — tt 
puzzle and needs more theoretical and experimental tests indeed. In this work, we consider 
the two possibilities and obtain the branching ratio of the decay X(3872) — ;■ ijj'^ with and 
without considering the mixing scenario respectively. 

In this work after the introduction we derive the form factors for the radiative decay 
of 2 ^ —7- 1 7 in the covariant light-front approach in section II. Then in section III we 
present our numerical results about the decay X(3872) — ;• J/ipl and X(3872) — )• ip'^y. The 
section IV is devoted to discussions and our conclusion. 

II. THE DECAY OF 2 + -^ 1 7 IN THE COVARIANT LIGHT-FRONT AP- 
PROACH 



In Refs. |36l . |37j the authors discussed how to calculate the transition matrix in the co- 
variant light-front approach. Following their strategy we formulate the matrix element for 
2 '" — >■ 1 7. The relative orbital angular momenta of 2^"^ and 1 are L = 2(d wave) and 



L = 0(s wave) respectively. 

Fist let us list the vertex functions for 2" 



45 



and 1" 



371 states as 



'i'H(^iD2)l5Kf,K^, 



iHv[y^ 



Wv 



(Pl - V2, 



MJ- 



where V represents the 1 state and K = ^2-^. The amplitude of ^D2 state decaying into 
1 via a photon emission is written as 



i 



pi 



-P2 



-V2 



FIG. 1: Feynnian diagrams depicting the radiative decay in the hght-front quark model. 
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^lut^u = Tr{75(-|^2 + m2)Yiu ^ ,^, ^ '' K/i + "^i)7/.(l^i + mi)}K^Kp, 



^afifiu 



Wv 



Tr{75(-^2 + "^2)7a«(-/2 + 'nT'2)[lu 



(Pl - P2)u 

Wv 



]i^i + nil)} KaK^ 
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Ni = p1 — m\ + ze, N[ = p\ — ml + ie, N2 = p^ — rn^ + ie and A^2 — v' 2 ~ ^^2 + ^^- The 
momentum pj is decomposed as [ vT ■ . Vi , Vi 1 ) in the light-front frame. One needs to integrate 
over pj~ by a contour integration [36|, l37|. Here we have set the relations among all relevant 



momenta a.s V = P + P', P = Pi + P2, P' = p'l + P2 Q = P — P', P2 = P2 for Fig i^)] Pi = P'l 
for Fig (b). Then it is easy to find s^^^^ = s^^^^. The integration contour is closed in the 
upper plane of the complex p^ for the first term in Eq.(2) and in the lower plane for the 
second term. Then the first term after the integration turns into 

where 



\^D2) 



^- (M'^-M'S)/M_L_,V.)for3., 
{M''-M'l)\ ^ J-. ct>\nD) ioi'D,, 

d) inD) = ^ — (hinb). 

In the formula, Wy and pi^, pi^pi^, Pi^PuPia ^^it^PiiyPiaPip in s;^^^^ must be replaced by 
the appropriate quantities as discussed in Ref.[37| to include the contributions of the zero 
modes, for example 

Wv -^ wv = Mq + nil + 7712 for ^5*1, 
Mq — {nil + "^2^"^ 



2Mq + mi + ni2 
. ^i-n , (^ P±Q±^ 



for ■'Di, 



with i°'^ and e' being identical to e"'^ and 5"^ for the maximally transverse polarized state 



[m 



:J) and s^^^j, changes into s^^^j,. More details about the derivation can be found in 
Ref.[37l|. After the replacement, s^/j^;^ is decomposed into 



with 



Sa^t^u = Fi {ep^p^VPq'^g^, + e^^p^pPq^^gp,) + {F^ + F^)e,^p^V'q'^q^qp 

+F^{epf,p^Vf'q'^qaqu + eaf.p^V^q'^qpq,), (4) 



4^AiW 



F2 = 2^^4(2) +22^2(2) -4^3(2)^ 



^ 4zA4(^) - 4^^4(2) ^ 8^^g(4) 
Wv 

Fa = . 

Wv 

where Ai^^\i = 1 ~ 4, j = 1 ~ 4) are defined in the appendix. 
We define the form factors as following 



/i 

h 
h 



167r3 



/ dxicfp± 



ee 



ecg 
ee„ 



'^ / dxi(fp± 
/ dxicfp± 



/ dxid'^p± 



Fi00' X1+X2, 

V^P^MoM'o a;iX2 ' 

F20(/)' a;i+a;2, 

^400' X1+X2, 



which will be numerically evaluated in next section. 
With these form factors the amplitude is obtained as 



(5) 



(6) 



^r\.n 



\h{ep^p^VPq''gau + e^^p^VPq'^gpu) + (/2 + h)eu^puV''q''qaqp 



+f\{.ep^p^V^q'^qo,qy + ea,^p^V^'q'^qpqy)]e''e 



u aP 



(7) 



III. NUMERICAL RESULT FOR X(3872) -^ J/V'7 AND ^(25)7 

With the formulas derived in section II we can calculate the decay rates for X(3872) — t- 
'0(15', 25')7 in the light-front quark model. In fact the main task is to evaluate these form 
factor /i, /2, /s and /4 at g^ = 0. 

We set TTic = 1.4 GeV following Ref. 37| at first and will discuss the dependence of our 



results on the choice in later part of the paper. Then we need to fix the parameter /3 in the 
wavefunction. Generally in LFQM, we can fix the /3 value by fitting the decay constant of a 
vector meson, thus in this work, we employ J /ip as the meson. The formula for calculating 
the decay constant of a vector meson was given in Refs. 36|, l37l |: 

<i){nS) 



h 






xMl 



dxi d'^p± 



?7li(?Tii — 7722) — Pj_ + 



2a;(l - x)Mo 

mi + 1712 



-pI 



Mo + mi + 1712 

where rrii = 7712 = rric and other notations are collected in the appendix. 

Using the data of J/ip — )■ e+e"!^?! we extract the decay constant of J/ip: fj/^ = 416 ± 5 
MeV by which we fix the model parameter P as 0.631 ± O.OOSGeV. The masses of X(3872), 
J/-0 and ■0(25') are chosen according to Ref.[47|. 
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TABLE I: the form factors /i, /2, /s, and f^ 



decay mode 



/i 



/2 



/3 



U 



X(3872) -^ J/'i/'T 

X(3872) ^ ^^(25)7 

X(3872) -^ ^{2S - lD)j 



-0.0146 ±0.0002 
-0.0157±0.0002 
-0.0594±0.0004 



0.146±0.003 
0.342±0.006 
0.330±0.005 



-0.0092±0.0001 
-0.0144±0.0004 
-0.0396±0.0007 



O.OlSOib 0.0001 
0.0352±0.0002 
0.0800±0.005 



The numerical values of the form factors for X(3872) — )■ J/ip'j are listed in Tab. I. Then 
we obtain r(X(3872) -^ J/i'l) = (3.54 ± 0.12) x lO^^GeV. The recent measurement sets 
only an upper bound on the total width of X(3872), so that we cannot determine accurate 
branching ratios for the radiative decays, but only a bound which is listed in table lllli 

Theoretically, Jia et.al. 2M studied the same transition within the pNRQCD framework 
where several phenomenological potentials were adopted, in their work the quantum num- 
ber of X(3872) was assigned as 2"+. They obtained r(X(3872) ^ J/^7) = (3.11 ~ 
4.78) X 10~^GeV corresponding to different potential models which is in accordance with 
our result. For comparing theoretical results with the data we present the experimental data 
and theoretical prediction on the branching ratios BTZ{X{3872) — ?> J/ip'~f) in table IIIII The 
theoretical prediction of BTZ{X{3872) — >■ J/ip'y) is slightly lower than the Babar's data but 
quite close to the Belle's data as long as we suppose X(3872) as a 2 ^ state. 

Now we turn to X(3872) — )■ ip''-/. We will deal with the cases where ijj' is regarded as 
a pure 25* charmonium ip{2S) 47| or a mixture of 2S and ID charmonium ip{2S — 1D)|46|, 



48 



49] respectively and obtain the corresponding numerical results. First we assume ip' 



to be a pure 2S state to calculate the branching ratio of the transition X(3872) — ;■ ilj''y. 
Here we need to use the wave function for the radial excited state. In Ref. 4lJ we notice 
that keeping the orthogonality among the nS states of heavy quarkonia and requiring the 
theoretical predictions on the decay constants of heavy quarkonia to be in agreement with 
that obtained from the data of leptonic decays, the wavefunction of the radial excited states 
nS {n >1) should be modified 4l|. In Eq.A(2) of the appendix we present the explicit 
form of the modified 25* wavefunction. With those wavefunction we obtain the form factors 
for X(3872) — !■ ilj{2S)'~f which are listed in Tab. I. Now we proceed to evaluate the decay 
rate of X(3872) -^ ^(2^)7. We obtain r(X(3872) -^ ^(2^)7) = (2.44 ± 0.10) x lO-^GeV. 
By the same assumption, the authors of Ref. 2J| obtained r(X(3872) — t- ilj(2S)'~f) = 1.7 ~ 
2.9 X lO-^GeV. 

The predicted branching ratio of B7l{X (3872) — )■ ^(2S')7) presented in table IIIII is much 
lower than the Babar's data, but does not contradict to the upper bound set by the Belle 
data. 



The authors of Ref. 2^ noticed that the mixing of ip{2S) and tp{lD) in ip' can enhance 
the rate of r(X(3872) — )■ ip'''y) remarkably. In this work, we also evaluate the impact of the 
mixing in ip' on r(X(3872) — ;■ ip''j) with the same mixing scheme and mixing angle 6 = 12° 



TABLE II: Theoretical predictions on the decay widths of the radiative decay of X(3872), which 
is regarded as a ^D2 charmonium 



r(X(3872) -^ J/^7) 

r(X(3872) -^ ^(2S)7) 

r(X(3872) -^ V(25 - lD)-f) 



(3.11 ~ 4.78) X lO-^GeV 

1.7 ~ 2.9 X lO^'^GeV 

(4.9 ~ 5.6) X lO-'^GeV 



our results 



(3.54 ±0.12) X lO-^GeV 
(2.44 ±0.10) X lO-^GeV 
(3.65 ±0.11) X lO^'^GeV 



TABLE III: Experimental data and theoretical predictions on the branching ratios of the radiative 
decay of X(3872), which is regarded as a ^D2 charmonium 



[29] 



> (1.35 ~ 2.08) X 10-3 

> (7.39 ~ 12.6) X 10-6 

> (2.13 ~ 2.43) X 10"^ 



our results" 



Bn{X{3872) -^ J/i;-f) 

B7^(X(3872) -^ ^(25)7) 

Bn{X{3872) -^ V(25 - lD)-f) 



> 8.75 X 10-3 

> 2.97 X 10-2 



> 5.56 X 10-3 
~ 1.01 X 10-2 



> 1.59 X 10-3 

> 10.6 X 10-6 

> 1.59 X 10-^ 



" the total width r(X(3872)) < 2.3 MeV[47l are used 



as given m 



24| . We calculate the form factors which are also presented in table I and obtain 
r(X(3872) -^ ip{2S - 1L>)7) = (3.65 ± 0.11) x IQ-^GeV. Our result is in accordance with 



that of Ref. 2J] which is within a range of (4.9 ~ 5.6) x IQ-^GeV. One can notice that 
even though the mixing of 2S and ID is taken into account the theoretical prediction on 
the branching ratio of r(X(3872) — )■ ^'7) is still lower than the Babar's data by two orders. 



Moreover, we obtain the ratios 



g7^(X(3872)^V)(25)7) 



0.069 and 



g7^(X(3872)^i/i(2S-lD)7) 
B7^(X(3872)-i>J/V'7) 



0.1 



B'R,(X(3872)-5>J/^7) 

which obviously contradict to the Babar's data. 

Now we turn to study the dependence of the theoretical predictions on rric for the three 
radiative decay modes. We let rric vary within a reasonable range from 1.3GeV to 1.6GeV 
and see how the widths of r(X(3872) -^ J/ip-f), r(X(3872) -^ ^(2^)7) and r(X(3872) -^ 
%1){2S — lD)'-f) change. Indeed, it is noted that when we do so, we need to refit the (3 
parameter for different rric values using Eq.([H]). The dependence of our results on rric is shown 
in Figj2l For example, when rric varies from 1.3GeV to 1.6GeV the theoretical prediction 
on r(X(3872) -^ ^{2S - lD)-f) changes from 4.21 x 10-^ to 2.82 x 10-^ It is about 33% 
change, but the order of magnitude remains unchanged. The main goal of the work is to 
look for a criterion to determine the spin-orbit-identity of X(3872). As aforementioned, the 
Babar datum is several orders larger than the prediction, thus the theoretical uncertainty 
about rric does not influence the qualitative conclusion. 
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FIG. 2: The dependence of r(X(3872) -^ J/ip-y) and r(X(3872) -^ ^'7) on m^ 

IV. DISCUSSION AND OUR CONCLUSION 

Determining the structures of the newly observed resonances is an interesting and difficult 
job, but the study is highly rewarding, because one can eventually identify if they are exotic 
states. Retrospecting the history of searching glueballs at the charm energy range, the 
results are indeed discouraging. The new discovery offers an opportunity to recognize the 
hadrons which are not existing in the basic representations of SU{3), but as exotic states 
or even mixtures of glueball and other states. As a matter of fact, the charm energy range 
may be an abundant mine for such new exotic resonances. However, before one can identify 
any of them as an exotic state, he has to thoroughly investigate if it is an excited state of 
regular charmonium no matter it is radially or orbitally excited. X(3872) is such an unknown 
resonance. Many authors consider X(3872) to be a molecular state or tetraquark, instead, 
it is also suggested that it may be a 2 '" charmonium. 

In this work we calculate the rates of the radiative decays of X(3872) in the LFQM. We 
ffist deduce the amplitude and corresponding form factors in the amplitude of 2 ^ — )■ 1 7 
with the LFQM. Then we suppose the X(3872) is a pure 2""*" charmonium and calculate 
the widths of X(3872) -^ J/ip-f and X(3872) -^ ip'j. We obtain r(X(3872) -^ J/ip-f) = 
(3.54 ±0.12) X lO-^GeV, r(X(3872) -^ '^{2S)-f) = (2.44 ±0.10) x lO-^GeV, r(X(3872) -^ 
iIj{2S - 1D)7) = (3.65 ± 0.11) x 10-^GeV. 

Our theoretical predictions on r(X(3872) — )■ ^'7) is obviously inconsistent with the data 
of the Babar collaboration, but that on i37^(X(3872) — )■ J/^7) is consistent with data within 
the error tolerance. 

As a matter of the fact, the transition matrix element eventually depends on an overlap- 



9 



ping integral of the wavef unctions of the initial and final hadrons. Thus the node structure 
of the wavefunctions is extremely important. Generally, the state with the principal quan- 
tum number n possesses n — 1 nodes, thus the overlapping integral between X(3872) which 
is supposed to be an orbital excited state with n = 1 does not possess a node, and nei- 
ther J/ip{lS), but il>{2S) does have a node. Therefore, the integrand in the overlapping 
integration of X(3872) and J/ip is always positive, instead, in the overlapping integration 
between wavefunctions of X(3872) and x(:(2S), the integrand is positive on the left side of the 
node of the 2S state and negative on the right side of the node. Thus a cancellation effect 
exists. Therefore by a common sense, the transition X(3872) — ;■ ip{2S) + 7 is suppressed 
and 5(X(3872) -^ i){2S) + 7) should be smaller than 5(X(3872) -^ J j^) + 7). When the 
25 — \D mixing for ?/;' is considered, the overlapping integration is decomposed into two 
integrations as 

/, = ^ fdxid^pi, ^ ^i - ( '^'^^' )0[cosg0^(25) - sin^0'(lD)] 

where (j)'{2S), (f)'{lD), are the radial parts of the wavefunctions for ijj{2S), ip{lD) which are 
the ingredients in ip', and X(3872), 6 is the mixing angle of 25" and ID states in ip'. There 
is a node in (j)'{2S), but not in (f)'{lD), so the first overlapping integration is suppressed as 
discussed before, but not for the second one. Therefore, the predicted branching ratio is 
enhanced in comparison with the pure 25" state assignment for ip'. 

The data of the Babar and Belle collaborations on X(3872) — >■ z/;' -|- 7 are obviously 
different, so that more accurate measurements are needed. 

Our conclusion is that if the branching ratio of X(3872) -^ ijj' + '-/ is indeed as large as the 
Babar group measured, the 2~~^ charmonium assignment of X(3872) should be ruled out. 
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Appendix A: Notations 

Here we list some variables appearing in the context. The incoming meson in Fig. [T]has 
the momentum P = Pi + P2 where pi and p2 are the momenta of the off-shell quark and 
antiquark and 

pf = XiP"*", P2 = X2P^-, 

Pl± = XiP± + p±, P2± = X2P± - P±, (Al) 
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with Xi and p± are internal variables and xi + X2 = 1. 
The variables Mo, Mq and A'^i are defined as 

^2 ^ pj + mj ^ pl + ml 
° xi 0:2' 

Mo = \/M2-(mi-m2)2, 



^/OON ./^X^/"^ /^P^ / 2'^P^+Pi\/ > Pz+Pix /AON 



.2 , „2 



with p^ = ^ - ^^, 5 = 1/1.82, 02 = 1.88684 and 62 = 1.54943. 
The Aij{i = 1 ~ 4, j = 1 ~ 4) are 

A, (1) - ^ 4 Ji) - 4, (1) - P^-^^ 4, (2) - _„2 _ (P^ • g^)' 

\2 ^1 



^2(2) = (Ai«)2, ^3(2) = ^1(2^2(2), ^4(2) = (A^^'^y 



9 ' 



Al(3) = Ai«Ai2, ^2^=^) = A2(^Mi(2), ^3^=^) = Ai«A2(2\ ^4^3) = A2^'^A2^'\ 

3 q 
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